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Abstract 

We introduce the beta generalized exponential distribution that includes the 
beta exponential and generalized exponential distributions as special cases. We 
provide a comprehensive mathematical treatment of this distribution. We derive 
the moment generating function and the rth moment thus generalizing some re- 
sults in the literature. Expressions for the density, moment generating function 
and rth moment of the order statistics also are obtained. We discuss estimation of 
the parameters by maximum likelihood and provide the information matrix. We 
observe in one application to real data set that this model is quite flexible and can 
be used quite effectively in analyzing positive data in place of the beta exponential 
and generalized exponential distributions. 

keywords: Beta exponential distribution, Information matrix. Generalized ex- 
ponential distribution. Maximum likelihood estimation. 

1 Introduction 



Gupta and Kundu [1] defined the cumulative distribution function (cdf ) of the generalized 
exponential (GE) distribution by 

^^.^(x) = (1 - e-"^)-, x>0. (1) 

The two parameters of the GE distribution represent the shape {a > 0) and the scale 
parameter (A > 0) like the gamma and Weibull distributions. The distribution ([T]) is also 
named the exponentiated exponential distribution. Clearly, the exponential distribution 
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is a particular case of the GE distribution when a = 1. Gupta and Kundu also investi- 
gated some of their properties. The exponentiated WeibuU (EW) distribution introduced 
by Mudholkar et al. [7] extends the GE distribution and was also studied by Mudholkar 
et al. |8j, Mudholkar and Hutson [9J and Nassar and Eissa [13j. Nadarajah and Kotz 
[12] introduced four more exponentiated type distributions: the exponentiated gamma, 
exponentiated WeibuU, exponentiated Gumbel and exponentiated Frechet distributions 
by generalizing the gamma, WeibuU, Gumbel and Frechet distributions in the same way 
that the GE distribution extends the exponential distribution. They also provide some 
mathematical properties for each exponentiated distribution. 

The GE density function varies significantly depending on the shape parameter a. 
Also, the hazard function is a non-decreasing function if a > 1, and it is a non- increasing 
function if a < 1. For a = 1, it is constant. The GE distribution has lots of prop- 
erties which are quite similar in nature to those of the gamma distribution but it has 
explicit expressions for the distribution and survival functions like a WeibuU distribution. 
The gamma, WeibuU and GE distributions extend the exponential distribution but in 
different ways. Therefore, it can be used as an alternative to the WeibuU and gamma 
distributions and in some situations it might work better in terms of fitting than the 
other two distributions although it can not be guaranteed. Moreover, it is well known 
that the gamma distribution has certain advantages compared to the WeibuU in terms of 
the faster convergence of the maximum likelihood estimates (MLEs). It is expected that 
the GE distribution also should enjoy those properties. 

Consider starting from the cdf G[x) of a random variable, Eugene et al. |2] defined a 
class of generalized distributions from it given by 

1 rG{x) 

Fi^) = -^r^ uj'^-'ii-ujf-'du, (2) 



B{a,b) 







where a > and 6 > are two additional parameters whose role is to introduce skewness 
and to vary taU weight and -B(a, b) = 10°"^^ {1 — ioY'^duj is the beta function. The cdf 
G{x) could be quite arbitrary and F is named the beta G distribution. Application of 
X = G"\V) to 1/ ~ B{a, b) yields X with cdf 

Eugene et al. [2] defined the beta normal (BN) distribution by taking G{x) in to be 
the cdf of the normal distribution and derived some first moments. General expressions for 
the moments of the BN distribution were derived by Gupta and Nadarajah Nadarajah 
and Kotz [lO] considered the beta Gumbel (BG) distribution by taking G{x) as the 
cdf of the Gumbel distribution and provided closed-form expressions for the moments, 
the asymptotic distribution of the extreme order statistics and discussed the maximum 
likelihood estimation procedure. Also, Nadarajah and Kotz [TT] worked with the beta 
exponential (BE) distribution and obtained the moment generating function, the first 
four cumulants, the asymptotic distribution of the extreme order statistics and discussed 
the maximum likelihood estimation. Some of the Nadarajah and Kotz's pi] results were 
generalized by Cordeiro et al. [1] who considered the beta WeibuU distribution. They 
derived the moment generating function, the moments and the information matrix. 

We can write by 

= JG(.)(a,6), (3) 
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where Iy{a,b) = B{a,b)^^ w"'^^{l — w)^^^dw denotes the incomplete beta function 
ratio, i.e., the cdf of the beta distribution with parameters a and b. For general a and b, 
we can express (E]) in terms of the well-known hypergeometric function defined by 
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i=o ^ ^ ' 

where (a)j = a{a + 1) . . . (a + i — 1) denotes the ascending factorial. We obtain 

a n[a, b) 

The properties of F{x) for any beta G distribution defined from a parent G{x) in ([2]) 
could, in principle, follow from the properties of the hypergeometric function which are 
well established in the literature; see, for example. Section 9.1 of Gradshteyn and Ryzhik 

The probability density function (pdf) corresponding to ([2]) can be put in the form 

fix) = ^^fi^r'i^ - Gix)}'-'gix), (4) 

where we noted that f{x) will be most tractable when the cdf G{x) and the pdf g{x) = 
dG{x)/dx have simple analytic expressions. Except for some special choices for G{x) in 
([2]) as is the case when G{x) is given by ([I]), it would appear that the pdf /(x) will be 
difficult to deal with. 

We now introduce the four parameter beta generalized exponential (BGE) distribution 
by taking G{x) in ^ to be the cdf ([T]). The cdf of the BGE distribution is then 

F{x) = — i- ' ^ u^-\l - cof-'dco, x>0, (5) 
B{a,b} Jo 

for a > 0, 6 > 0, A > and a > 0. The pdf and the hazard rate function of the new 
distribution are, respectively, 

fix) = j^^e-'%l-e-'T^-'{l-{l-e~'T}'-\ ^ > 0, (6) 

and 

aAe-^^(l - e-^^)°'^-i |l - (1 - e"^^)"}^"^ 

= p. ,x J ^ , x>0. (7) 

B[a, b) ii_(^i_^-\x-^a[a, b) 

The density function ([6]) does not involve any complicated function. If X is a random 
variable with pdf ([6]), we write X ~ BGE{a,b, X,a). The BGE distribution generalizes 
some well-known distributions in the literature. The GE distribution is a special case 
for the choice a = b = 1. If in addition a = 1, we obtain the exponential distribution 
with parameter A. The BE distribution is obtained from ([5]) with a = 1. It is evident 
that (E]) is much more flexible than the GE and BE distributions. Plots of the density 
(E]) and failure rate ([7]) for some special values of a, b, A and a are given in Figures 
[1] and [21 respectively. A characteristic of the BGE distribution is that its failure rate 
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BGE(a,b,1,1) 



BGE(a,b,1,2) 




Figure 1: Plots of the density ([6]) for some values of the parameters. 



can be bathtub shaped, monotonically increasing or decreasing and upside-down bathtub 
depending basically on the values of its parameters. 

We now introduce the double generalized exponential (DGE) distribution for the 
choice a = 1 in ([5]). The cdf of the DGE distribution is 

= {l-(l-e-^")"}^ x>0. (8) 



The rest of the paper is organized as follows. In Section 2, we derive some expressions 
for the cdf of the BGE distribution and for the pdfs of the order statistics. In Sections 3 
and 4 we provide expressions for the moment generating function and for the rth moment, 
respectively. Furthermore, in these sections, we derive corresponding expressions for the 
order statistics. In Section 5, we discuss maximum likelihood estimation, inference on 
the parameters, and calculate the information matrix. One application to real data set 
is presented in Section 6. Some conclusions are drawn in Section 7. 
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Figure 2: Plots of the hazard ([7]) for some values of the parameters. 
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2 Distribution function and order statistics 



We provide two simple formulae for the cdf of the BGE distribution depending if the 
parameter 6 > is real non-integer or integer. First, if |^| < 1 and 6 > is real non- 
integer, we have 

Using the expansion (IHl) in ([5]), the cdf of the BGE distribution when 6 > is real 
non-integer follows 

Fix) = I u^-Hl - ujf-^dcu 



B{a,b) 



B{a,b)j^J{b-j)j\ 







and then 



^ ^ r(a) ^ Tib - ma + J) ■ ^^'^ 

Equation (fTUl) reveals the property that the cdf of the BGE distribution can be expressed 
as an infinite weighted sum of cdfs of GE distributions 

T{a + b) ^ (-l)^'GA,a(a+i)(x) 

^""^ r(a) 2-.r{b-ma + j)- 

Using the binomial expansion in ([5]), we have when 6 > is integer 

Fix) = — 1- E f ^ " ^) ^(1 - e-'T^'^^^\ (11) 

and, again, the same property of equation ffTUl) holds but now the sum is finite. Expres- 
sions ( fTOj) and ( fTTl) are the main results of this section. The cdf of the BE distribution 
follows with a = 1 from (fTOl) and ( ITTi) depending if 6 > is real non-integer and integer, 
respectively. The cdf of the DGE distribution follows from the above expressions when 
a = 1. The cdf ([T]) of the GE distribution comes from (llip when a = b = 1. If, in 
addition, a = 1, (fTTl) yields the cdf = 1 — e~'^^ of the exponential distribution. 
It can be seen in the Wolfram Functions that for integer a 



(i-!/)'^r(6 + j) 



m ^ J' 



and for integer b 



6-1 



IyM = --}_^^^il~yy. 



http : //functions .wolf reun. com/ 
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Therefore, for a > integer, 

and for 6 > integer we have an alternative form for (ITTi) given by 

r(a) ^ j! 

The density function (I6l) can be expressed in the mixture form in terms of cdfs of the 
GE distributions 

for b real non-integer and 
for 6 integer. 

Simulation of the BGE distribution is easy: if V is a random variable with a beta 
distribution with parameters a and b, then X = — log(l — V^^")/\ follows the BGE 
distribution with parameters a, b, A and a. 

We now give the density of the ith order statistic Xi-n, fi:n{x) say, in a random sample 
of size n from the BGE distribution. It is well known that 

hn{x) = ^ f{x)F^-\x) {1 - F{x)r-' , 

B[i, n — t + 1) 

for i = 1, . . . ,n. Using ([3]) and (Ej) we can express fr.n{x) in terms of the hypergeometric 
functions by 

. , aXe-^^G^,a^.l{x)Gx,aa(^-l){x){l - G x,^{x)}'^''-'+''>-' 
^""^"^^ B{i,n-i + l)B{a,b)^a'-^b^-' ^ 

2Fi(a, 1 - 6, a + 1; Gx,^{x)y-' 2^1 1 - a, 6 + 1; 1 - GA,.(a;))"-\ 

We now derive two alternative expressions for the densities of the order statistics using 
the expansion (Yll^i ^i)^ — Yl^Li mk}=o • • • ^ ^ positive integer. With this 

expansion and from ffTUj) — ffTTj) . we can show for 6 > real non- integer and integer that 

n—i 00 00 
k=0 mi=0 mfc+i_i=0 
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and 

n—i 6—1 b—1 

/-(^) = EE--- E (13) 

fc=0 mi=0 m(;+i_i=0 

respectively, where from now on fk,i{x) represents the density of a random variable Xk^i 
following a BGE{a{a{i + 1) + rrij} ,b, X, a) distribution, and the functions d^^l 

(2) 

and 5^ ■ required for the above expressions are 



and 



5(a,6)^+*S(i,n + 1) _^ 

The sums in (fT^ and (fT5I) extends over all (fc + 2)-tuples {k,mi, . . . ,mfc+j„i) of non- 
negative integers and are easily implementable on a computer. 

3 Moment generating function 

The moment generating function (mgf) of the BGE distribution is given by 



aX 



oo 



M{t) = ^7'' . / e*^e-^^(l - e-^^)"'^-^! - (1 - e-^^)"}^-Ma;. (14) 
B{a, b) Jq 

Using the expansion when 6 > is real non-integer, ffT^ reduces to 
Changing to the variable u = e~^^, we obtain 

«r(6) ^ (-1)^- r 

i=o 



The above expression shows that the mgf of the BGE distribution exists if t < A. As- 
suming that t < A, we have 

Analogously, if 6 > is integer, assuming t < A, we use the binomial expansion in ffT^ 
to obtain 

M{t) = ^ f ~ ^) (-l)'^(l - VA, «(a + j))- (16) 



If we take a = 6 = 1 in (fT6l) . the above mgf reduces to 

M{t) = aB{l-t/\,a), 

which agrees with Gupta and Kundu's [5] equation (2.3). 

From ( fT5i) and ( fT6l) with a = 1 we have the mgf of the BE distribution 

B{a, b) 

which agrees with Nadarajah and Kotz's [IT] equation (3.1). The last equahty is given 
with more details in the Appendix A. The expressions for the mgf of the DGE distribution 
are obtained with a = 1 in ( |T5l) and ( |T6l) for 6 > real non-integer and 6 > integer, 
respectively. 

We can easily obtain expressions for the mgf of the order statistics from (fT^ and (fT^ . 
The mgf of for 6 > real non-integer is 

n—i oo 
fc=0 mi=0 

and for 6 > integer 

n—i 6—1 
k=0 mi=0 

where Mk^i{t) is the mgf of the BGE{a{a{i + 1) + Xl^i"^ "^j}-, ^5 -^^ oi) distribution. 



mfc+j_i=0 



mfe+,_i=0 



4 Moments 

The rth moment of the BGE distribution can be obtained from E{X^) = dJ" Mx{t) / dt'^\t=Q. 
Hence, if 6 > is real non-integer, we have from 0151) 



aT{b) 

XrB{a,b)j^^T{b-j)j\ dp- 



p=i 



and if 6 > is integer, we obtain from (|T6 



aT{b) 
\'B{a,b) 



6-1 



6-1 



d'-B{p, a{a + j)) 



dp- 



(19) 



Equations ( fT8|) and (fT9l) are the main results of this section and generalize the moments 
of the GE and BE distributions derived by Gupta and Kundu [5] and Nadarajah and 
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Kotz [TT], respectively. The first four moments of tlie BGE distribution for 6 > real 
non-integer are 

, T{a + h) ^ {-iy{a + 3)-^ , r(a + &) ^ (-l)^(a + jQ-^ 

Ar(a) ^ V{h-j)j\ ^2-pr(a)^ V{h-j)j\ ""^^ 

^ _ T[a + h) ^ {-iy+\a + 3)-^ , _ r(a + &) ^ (-l)^(a + 

A3r(a) ^ r(6-j)j! A4r(a) ^ r(6-j)j! 

and for 6 > integer are 

Ar(a) ^1, J Ja + /- A2r(a) ^1, j Ja + /- 

, _r(a + 6)^ , r(a + &) ^/6-l\(-l)^- 

A3r(a) J ) a + j Xna) J Ja + /- 

where the quantities Cj, dj, ej and fj are given in the Appendix B. 

Graphical representation of skewness and kurtosis for some choices of parameter b as 
function of parameter a, and for some choices of parameter a as function of parameter b, 
for fixed A = 1 and a = 1, are given in Figures [3] and HI respectively. These plots show 
that the skewness and kurtosis curves increase (decrease) with a (6) for fixed b (a). 

The rth moment of for 6 > real non-integer is 



n—i oo 



Li) 



k=0 7711=0 mk+i-i=0 

and for 6 > integer 

n—i b—1 b—1 
k=0 mi=0 ms;_|_i_i=0 

where the moments E{Xl^) come from the general expansions ffT^ and ffT^ for the mo- 
ments of the BGE distribution with parameters a{a{i + 1) + Yl'j=i~^ "^j}) ^) ^ ^ind a. 

The Shannon entropy of a random variable X is a measure of the uncertainty and 
is defined by log /(X)}, where f{x) is the fdp of the X. For a random variable X 
with a BGE distribution, we obtain 

E{-log/(X)} = -log(aA)+logS(a,6)+A/i; + (l/a-a){^(a)-^(a + 6)} 
-(6-l){^(6)-^(a + 6)}, 

where fi'i is the mean given before. 
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Figure 3: Skcwness of the BGE distribution as function of a for some values of b and as 
function of b for some values of a. 




Figure 4: Kurtosis of the BGE distribution as function of a for some values of b and as 
function of b for some values of a. 
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5 Estimation and Inference 



We assume that Y follows the BGE distribution and let 9 = (a, b, A, ce)^ be the parameter 
vector. The log-hkehhood for a single observation y of F is 

£ = e(a, b, A, a) = log a + log A - log B{a, b) - Xy + (aa - 1) log(l - e~^^) 

+(6-i)iog{i-(i-e-^T}, y>o. 

The components of the unit score vector U = (|^, ||, ||, are 



d£ 

- = alog(l 



e-^y)-^(a) + ^{a + b), 



m 

db 

de 

dl 
da 



log{l - (1 - e-^^)"} - V'(fc) + i^{a + 6), 



1 ^ ^ (aa - l)ye-^2/ ^{b - l)e-^y{l - e'^y^-^y 



1 - e-^y 1 - (1 - e-^yy 

{b - 1)(1 - e-^y)"log{l - e-^y) 



+ alog(l -e"^^') - 



a ' 1 - (1 - e-^y)"' 

The expected value of the score vector vanishes and then 

ip{a) — ip{a + b) 



E{\og{l - e-'"^)} 
E[\og{l - (1 - e-^^)"}] 



a 



E 



(1 



-\Y\a 



)"log(l -e-^^) 



1 - (1 - e-^^)" 



a{V>(a) --0(a + &)} + 1 
q;(6-1) 



For a random sample y = (yi, . . . of size n from Y , the total log- likelihood is = 
£„(a, 6, A,q;) = where ig the log-likelihood for the ith observation [i — 

The total score function is C/„ = S^Li ^^'^ where C/W has the form given 
before for i = 1, . . . ,n. The MLE ^ of ^ is obtained numerically from the nonlinear 
equations ?7„ = 0. For interval estimation and tests of hypotheses on the parameters in 
6 we obtain the 4x4 unit information matrix 



K = K{9) = 













K'b,b 




^b,a 




K.b,X 








t^b,a 
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where the corresponding elements are given by 

cx 

l^a,a = ^'{a) - ■ip'ia + b), Ka,b = {a + b), Ka,X = ^7o,i,i,i,o 

V>(a + 6) -jjja) a 

K'a,a = , K'b,b = W (o) - (a + o), K^^x = --Ti^l,l,l,0, 

CX A 

«A,A = ^[1 + («« - l)(ro,2,2,2,0 + ^0,1,1,2,0) + Oi{b - 1){q;T2,2,2,2,0 + 

T 11 ^ _ a{V'(a) --0(q + &)} + ! 
{a - ij J 1,2,2,2,0 - 1,1,1,2,0 jj, H,a ' 

«A,a = ^{aTo,i,i,i,o - (6 - l)(7l,i,i,i,o + ^2,1,1,1,1 + 
«a,a = ^{1 + l)(r2,o,0,0,2 + ri,o,o,0,2)}. 

Here, we have defined the following expectation 

where V ~ Beta{a,b) and i,j,k,l,m G {0,1,2}. The total information matrix is then 
Kn = Kn{e) = nK{e). 

Under conditions that are fulfilled for parameters in the interior of the parameter 
space but not on the boundary, the asymptotic distribution of 

^Jn{e-e) is N4{o,K{e)~^). 

The asymptotic multivariate normal N4^{0, Kn{0)~^) distribution of 9 can be used to 
construct approximate confidence intervals and confidence regions for the parameters 
and for the hazard and survival functions. The asymptotic normality is also useful for 
testing goodness of fit of the BGE distribution and for comparing this distribution with 
some of its special sub-models using one of the three well-known asymptotically equivalent 
test statistics - namely, the likelihood ratio (LR) statistic, Rao score (Sr) and Wald (W) 
statistics. 

An asymptotic confidence interval with significance level 7 for each parameter 9i is 
given by 

ACI{9i, 100(1 - 7)%) = {9i - z^i2^f^\ 9i + z^/2\f^^), 

where is the zth diagonal element of Kn(9)~^ for i — 1,2,3,4 and is the quantile 
1 — 7/2 of the standard normal distribution. 

Further, we can compute the maximum values of the unrestricted and restricted log- 
hkelihoods to construct the LR statistics for testing some sub-models of the BGE dis- 
tribution. For example, we may use the LR statistic to check if the fit using the BGE 
distribution is statistically "superior" to a fit using the GE distribution for a given data 
set. In any case, considering the partition 9 = {9i,92Y , tests of hypotheses of the 
type Hq : 9i = 9^^^ versus Ha '■ 9i ^ ^1°'* can be performed by using any of the above 
three asymptotically statistics. The LR statistic for testing the null hypothesis Hq is 
w — 2{i{9) — i{9)}, where 9 and 9 are the MLEs of 9 under Ha and Hq, respectively. 

Under the null hypothesis, w Xqj where q is the dimension of the vector 9i of interest. 
The LR test rejects Hq ii w > ^j, where denotes the upper 1007% point of the xl 
distribution. 
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6 Application 



In this section we fit BGE model to one real data set. The data set is obtained from 
Smith and Naylor [14J. The data are the strengths of 1.5 cm glass fibres, measured at 
the National Physical Laboratory, England. Unfortunately, the units of measurement are 
not given in the paper. The data set is: 0.55, 0.93, 1.25, 1.36, 1.49, 1.52, 1.58, 1.61, 1.64, 
1.68, 1.73, 1.81, 2 ,0.74, 1.04, 1.27, 1.39, 1.49, 1.53, 1.59, 1.61, 1.66, 1.68, 1.76, 1.82, 2.01, 
0.77, 1.11, 1.28, 1.42, 1.5, 1.54, 1.6, 1.62, 1.66, 1.69, 1.76, 1.84, 2.24, 0.81, 1.13, 1.29, 1.48, 
1.5, 1.55, 1.61, 1.62, 1.66, 1.7, 1.77, 1.84, 0.84, 1.24, 1.3, 1.48, 1.51, 1.55, 1.61, 1.63, 1.67, 
1.7, 1.78, 1.89. 

The MLEs of the parameters and the maximized log-likelihood for the BGE distribu- 
tion are 

a = 0.4125, b = 93.4655, A = 0.92271, a = 22.6124, isGE = -15.5995, 

whereas for the BE distribution 

d = 17.7786, b = 22.7222, A = 0.3898, Ibe = -24.1270 

and for the GE distribution 

A = 2.6105, a = 31.3032, Ige = -31.3834. 

The LR statistics to test the hypotheses Hq : BE x Ha '■ BGE and Hq : GE x Ha '■ 
BGE are 17.0550 (p-value=3.63 x lO^^) and 31.5678 (p-value=1.39 x 10"^), respectively. 
Therefore, we reject the null hypothesis in both cases in favor of the BGE distribution 
at the significance level of 5%. The plots of the estimated densities of the BGE, BE and 
GE distributions fitted to the data set given in Figure [5] show that the BGE distribution 
gives a better fit than the other two sub-models. 

7 Conclusions 

We proposed the beta generalized exponential (BGE) distribution which generalizes the 
beta exponential distribution discussed by Nadarajah and Kotz [Tl] and the generalized 
exponential (also named exponentiated exponential) distribution introduced by Gupta 
and Kundu 0]. We provide a mathematical treatment of this distribution including 
the densities of the order statistics. We derive the moment generating function and 
provide infinite sums for the moments of the new distribution and of its order statistics. 
We discuss maximum likelihood estimation and obtain the information matrix. One 
application of the BGE distribution are given to show that this distribution could give 
better fit than other sub-models discussed in the literature. We hope this generalization 
may attract wider applications in reliability and biology. 
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Appendix 

We now give some details of ffTTj) . It follows for t < A that 



»i 

B{h-t/\,a) = B{a,h-t/\)= I co^-^l - cuy-'^^duj 



3=0 



oo 



3=0 

oo 



usin 




(-i; 


^T(6) 


r(6- 




(-i; 


^T(6) 


r(6- 




(-i; 





Similarly for 6 > integer and using the binomial expansion in (1 — cj)*"^, we have 
Bib-t/X,a) = J2{ ■ lB{l-t/X,a + j). 

3=0 ^ ^ ^ 

Hence, it is easy to see that ( fTTll holds. 
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Some useful quantities are defined by 

Cj = V(a(a + j) + l)-^(l), rf, = c| + ^'(l)-^'(a(a + j) + l), 

= -c,[c2 + 3{V''(l)-V^'(a(a + j) + l)}]+^"(l)-^"(«(a + j) + l), 
/. = {c2 + vy(l)-^'(«(a + j) + l)}[s' + 3{^'(l)-^(a(a + j) + l)}] + 
2c|{V^'(l) - ij'{a{a + 3) + 1)} - 4c,{<(l) - <(«(a + j) + 1)}. 
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